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The likelihood that an undercooled liquid vitrifies or crystallizes depends on the cooling rate IZ. 

The critical cooling rate TZc, below which the liquid crystallizes upon cooling, characterizes the 
glass-forming ability (GFA) of the system. While pure metals are typically poor glass formers with 
IZc > K/s, specific multi-component alloys can form bulk metallic glasses (BMGs) even at 
cooling rates below 77. ~ 1 K/s. Conventional wisdom asserts that metal alloys with three or more 
components are better glass formers (with smaller IZc) than binary alloys. However, there is cur¬ 
rently no theoretical framework that provides quantitative predictions for IZc for multi-component 
alloys. In this manuscript, we perform simulations of ternary hard-sphere systems, which have been 
shown to be accurate models for the glass-forming ability of BMGs, to understand the roles of geo¬ 
metric frustration and demixing in determining IZc- Specifically, we compress ternary hard sphere 
mixtures into jammed packings and measure the critical compression rate, below which the system 
crystallizes, as a function of the diameter ratios cts/cta and ctc/ua and number fractions xa, xb, 
and xc- We find two distinct regimes for the GFA in parameter space for ternary hard spheres. 

When the diameter ratios are close to 1, such that the largest (A) and smallest (C) species are 
well-mixed, the GFA of ternary systems is no better than that of the optimal binary glass former. 

However, when acjaA / 0.8 is below the demixing threshold for binary systems, adding a third 
component B with ac < ctb < ua increases the GFA of the system by preventing demixing of A 
and C. Analysis of the available data from experimental studies indicates that most ternary BMGs 
are below the binary demixing threshold with acjoA < 0.8. 

PACS numbers: 64.70.pe,64.70.Q-,61.43.Fs,61.66.Dk,61.43.Dq 


1. INTRODUCTION 

When atomic and molecular liquids are cooled suffi¬ 
ciently rapidly (i.e. above the critical cooling rate TZc), 
they bypass crystallization and become trapped in dis¬ 
ordered glassy configurations [l|. Avoiding crystalliza¬ 
tion in pure metals is very challenging and has only been 
achieved in experiments recently • On the other hand, 
multi-component liquid alloys can form bulk metallic 
glasses (BMGs) that possess centimeter or greater casting 
thicknesses and critical cooling rates TZc < 1 K/s 
BMGs have shown great promise as structural materials 
because they are amorphous with few defects andpossess 
higher processability than crystalline metals @, [3- 

The conventional wisdom in the BMG research com¬ 
munity is that BMGs should contain three or more 
atomic species Q with atomic size differences above 12% 
{i.e. the ratio of the diameters of the smallest to the 
largest species should be < 0.89) Intuitively, more 
atomic components with different sizes introduces geo¬ 
metric frustration or “confusion,” which delays crystal¬ 
lization BEini. Also, it has been suggested that a 
mixture of multiple atomic species leads to dense pack¬ 


ing in the liquid state and thus enhanced stability of the 
glass [HI- The minimum critical cooling rate observed 
for binary BMGs is TZc ^ 10^ K/s, while it decreases 
to 10“^ and 10“^ K/s for ternary and quaternary sys¬ 
tems, respectively. Alloys with similarly sized atomic 
constituents can only be cast into glassy thin films. (See 
Table [Jin Appendix lAl) Even with these empirical rules, 
there is an enormous parameter space of potential BMGs 
and we lack a complete theoretical framework that would 
enable the prediction of TZc for each alloy in the design 
space. 

A number of recent studies of hard-sphere mixtures 
have shown that dense atomic packing plays an impor¬ 
tant role in determining the glass-forming ability (GFA) 
of metallic alloys [12 - ll7l| . In particular, binary metal- 
metal {i.e. transition metal-transition metal) BMGs, 
such as Cu-Zr, Cu-Hf, and Ga-Al, possess atomic size 
ratios a = obIoa and compositions Xb that occur in 
the region of parameter space with the smallest TZc for 
binary hard spheres [l8l |. 

Gan hard sphere models accurately capture the depen¬ 
dence of the GFA on the atomic size ratios and compo¬ 
sitions for multi-component alloys? In this manuscript. 
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FIG. 1: (Color online) (top) Comparison of the volumes of 
tetrahedral cells formed by the centers of face-centered cu¬ 
bic packed atoms with (left) two versus (right) three different 
sizes {i.e. either gcIoa = 0.8 with xb = 0 or obIga = 0.9, 
gcIga = 0.8, and xb > 0). The tetrahedral volume for four 
same-sized atoms is 0.118 g\ (thin black line), while it is 
0.084 ga for two large and two small atoms (thick blue solid 
line) and 0.092 ga for one small, one intermediate, and two 
large sized atoms (red dashed line). Thus, the distortion of 
the tetrahedral cell is smaller for the ternary system, (bot¬ 
tom) When atoms of an intermediate size (gbIga = 0.75) are 
added to a binary system (with diameter ratio gc(ga = 0.5), 
it becomes more uniformly mixed and less ordered, i.e. with 
global bond orientational order parameter Qq = 0.02 (right) 
compared to 0.15 (left). 


we study the glass-forming ability of ternary hard-sphere 
mixtures. We find two key results: (1) When the sizes of 
the three components are comparable, ternary systems 
behave similar to binary systems, and the GFA cannot 
be larger than that of a binary system consisting of the 
largest and smallest components. In this case, the pack¬ 
ing efficiency of the ternary system is close to that of the 
binary systems. (See Fig. [T] (top).) (2) When the di¬ 
ameter ratio of the smallest to the largest component is 
beyond the demixing limit (a ;< 0.8 [l8j|). adding a third 
component with an intermediate size can increase the 
GFA by preventing demixing. In this scenario, the pack¬ 
ing fraction of the ternary system is significantly higher 
than the demixed binary system. (See Fig. [T] (bottom).) 
This demixing mechanism has also been found in studies 
of segregation of granular media and other particulate 
solids |1^ . 


2. METHODS 

We performed event-driven molecular dynamics sim¬ 
ulations oi N = 500 ternary hard spheres with diam¬ 
eters (JA > <^B > crc, number fractions xa = Na/N, 
Xb = Nb/N, and xc = Nc/N, and the same mass m. 
We compressed systems initially prepared in liquid states 


at packing fraction (j) = 0.25 so that they exponentially 
approach static jammed packings at ^ as a func¬ 

tion of time. We terminate each compression run when 
(^j — 0)/</>j = 10“^. We vary the compression rate R 


over 5 orders of magnitude m- Note that R is given in 
units of y/kBT/ma\ and R = 1 corresponds to a cool¬ 
ing rate TZ ^ 10^^ K/s for metal alloys [l^. The crys¬ 
tal structures that compete with glass formation possess 
face-centered cubic (FCC)-like order, and thus we char¬ 
acterize the positional order of the packings using the 
global bond orientational order parameter Qq aver¬ 
aged over 96 independent compression runs. The critical 
compression rate Rc is determined by the intersection 
of the mean and median Qq as a function of R. (See 
Appendix [BJ) To explore the glass-forming ability dia¬ 
gram for ternary systems, we studied more than 20 com¬ 
positions and 10 pairs of atomic size ratios ob/^a and 
A- Additional details of the simulation methods can 
be found in Ref. fisl. 


3. RESULTS 

In our previous studies of binary hard-sphere mixtures 
with diameter ratio a = gb I ga and small particle num¬ 
ber fraction XB^ we found that the critical compression 
rate Rc decreases exponentially, Rc ^ exp[C(a;s)(l — 
q:)^], where C is a composition-dependent constant, for 
a > 0.8 above the demixing limit. (See Fig. [^) In con¬ 
trast, for a < 0.8 the large and small particles in binary 
systems can demix, which then induces crystallization. 
Thus, the glass-forming ability for binary hard-sphere 
systems first increases with decreasing a, but then be¬ 
gins to decrease for a < 0.8 (l8l| . 

We first focus on ternary hard-sphere systems with 
weak size disparities. In Fig. [3] (top), we plot the critical 
compression rate Rc as a function of the number fractions 
of the three components xa, xb and xq at fixed diameter 
ratios gbIga = 0.95 and gcIga = 0.9. We find that the 
best glass-forming regime (i.e. the region with the small¬ 
est Rc) occurs on the binary composition line AC near 
Xc = 1 —~ 0.6 and xb = 0. Adding the third compo¬ 
nent B with an intermediate size gq < gb < ga causes a 
decrease in the glass-forming ability (or increase in Rc). 
Since the competing crystal for bidisperse hard-sphere 
systems in this diameter ratio regime is a deformed FGG 
crystal, adding a third component with an intermedi¬ 
ate size to a binary system reduces the lattice distortion 
(Fig. [D (top)) and the glass-forming ability. These re¬ 
sults are consistent with experimental observations for 
bulk metallic glasses, which are summarized in Table U in 
Appendix]^ There are no observed ternary bulk metallic 
glasses with weak size polydispersity, i.e. the diameter 
ratio of the smallest to the largest component satisfies 
a > 0.8. 

We now consider ternary hard-sphere systems with 
a diameter ratio disparity that is beyond the demixing 
threshold, i.e. gcIga < 0.8. In Fig. [3] (middle), we 
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FIG. 2: (Color online) In binary hard-sphere mixtures with 
diameter ratio a > ac{xB) and fixed number fraction xb of 
small particles, log^Q drops linearly with (1 — a)® (dotted 
line). For a < Oc, the large and small particles demix and 
Rc begins to increase with decreasing a. The composition- 
dependent threshold is Oc « 0.8 for xb = 0.2. The inset 
shows snapshots of configurations at R ^ Rc for xb = 0.2 
and a = 0.9, 0.7, and 0.5 from left to right, which illustrates 
increasing demixing as a decreases. 


plot the critical compression rate Rc as a function of the 
compositions xa, xb, and xc for ternary systems with 
diameter ratios (JbI<Ja = 0.95 and adoA = 0.5. For 
this system the smallest value of Rc does not occur for 
Xb = 0. In contrast, we show that Rc decreases with 
increasing xb- (See Fig. [3] (bottom).) 

We can also measure the glass-forming ability at fixed 
composition and vary the diameters of one of the parti¬ 
cles. In Fig. m we fix the compositions xa = xb = xc = 
1/3 and diameters a a and ac of two components and 
measure Rc as a function of the diameter ratio ubI<xa- 
Note that, when ob = oa the ternary systems 

reduce to binary systems with xc = 1/3 (2/3). In ex¬ 
perimentally observed ternary BMGs, when the diame¬ 
ters of two of the three components are similar, for in¬ 
stance CuNi and AlTi, the ternary glass-forming ability 
diagram is symmetric and equivalent to that of the cor¬ 
responding binary system [22| . We first focus on ternary 
systems with ac I <xa = 0.9, which does not lead to demix¬ 
ing. When acl<JA < <xbI<xa < !> Rc has a maximum at 
<xbI<xa < 1 and the ternary systems are worse glass for¬ 
mers than binary systems with ob = <Jc- These ternary 
systems show enhanced glass-forming ability above that 
for binary systems only when gb^a ^ 0.9. (See Fig. ID) 

In Fig. |4l we also consider fixed diameter ratio 
(tedA = 0.5 for which the two components tend to 
demix. In this case, Rc does not possess a maximum 
at aBl<JA < 1) and thus these ternary systems can pos¬ 
sess enhanced glass-forming ability compared to corre¬ 
sponding binary systems. The introduction of the third 
component with an intermediate size gb prevents demix¬ 
ing. As shown in the insets of Fig. 01 binary systems with 
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FIG. 3: (Color online) (top) Gritical compression rate Rc as 
a function of the compositions xa, xb, and xc in ternary 
hard-sphere systems with diameter ratios gbIga = 0.95 and 
GcjGA = 0.9. The minimum Rc occurs on the edge AC, i.e. 
binary systems with xs = 0. Moving perpendicular to AC on 
the diagram causes increases in Rc. (middle) Same as the top 
plot, but with diameter ratios gbIga = 0.95 and gcIga = 
0.5. The minimum Rc no longer occurs for binary systems 
on the edge AC. Note that the contour plots of logj^g Rc are 
interpolated from ~ 20 simulation runs (triangles) and given 
on a color scale that decreases from dark to light, (bottom) 
Rc as a function of xb at fixed xc = 0.4 (dashed line in top 
and middle panels) for the diameter ratios studied in the top 
(squares) and middle (circles) panels. 
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FIG. 4: (Color online) The critical compression rate Rc as 
a function of the diameter ratio gbIo-a at fixed composition 
XA=XB=xc = ^l‘i and acl^A = 0.9 (squares) or ac jcFA = 
0.5 (circles). The solid lines are polynomial fits to the data 
to show qualitative trends. Ternary systems reduce to binary 
systems when gb = ga or gc (solid symbols). In the inset, 
we show configurations obtained at a slow compression rate 
R = 10 “® for gbIga = 0.5 (left) and 0.75 (right). Large, 
intermediate, and small particles are shaded from dark to 
light. 


cts/cta = 0.5 demix and crystallize (left), while ternary 
systems with ub!( x a = 0.75 remain well-mixed and amor¬ 
phous (right). Although the large particles A exclude the 
small ones C, A particles mix with B particles and B 
particles mix with C particles, which leads to effective 
mixing of A and C particles. 

Because packing efficiency and vibrational entropy 
determine the stability of crystals in hard-sphere sys¬ 
tems Hi , one can correlate the packing fraction at 

jamming ij)j with the critical compression rate Rc as 
demonstrated in binary systems [l^. We study three 
relevant packing fractions: obtained in the limit 

R ^ oo [ 2 ^ [ 2 ^, (j)j for amorphous packings obtained 
a.t R ^ Rc, and (/>} for partially crystalline packings ob¬ 
tained aX R Rc- We find that the packing fraction 
of single-crystal FCC packings is not strongly correlated 
with Rc- Instead, partially crystalline systems compete 
with glass formation (Tsj l. 

As shown in Fig. [S] for ternary systems with diam¬ 
eter ratio pairs gbI^a = 0.95 and uc/ga = 0.9 and 
xbI<xa = 0.9 and gc/ga = 0.88 that do not demix, the 
relations between Rc and packing fraction 4>j at jamming 
follow the trends for binary systems, i.e. as (j)°j and (jjj 
approach each other, Rc -A 0. In addition, the pack¬ 
ing fraction in these ternary systems is not larger than 
in binary systems, contrary to the intuition that ternary 
systems are always denser than binar y sy stems and thus 
possess higher glass-forming abilities [Tl|. These results 
show that ternary systems with weak diameter ratio dis¬ 
parities can be described effectively as binary systems 
since the additional intermediate-sized particles only de- 



FIG. 5: (Golor online) (a) For each binary and ternary system, 
we plot the corresponding critical compression rate Rc and 
three definitions of the packing fraction at jamming: 
obtained in the i? —>■ 00 limit, (j)j for amorphous packings 
obtained at R Rc, and (j^j for partially crystalline packings 
obtained at A ~ Ac with < (j)°j < (jXj. The solid lines 

give polynomial fits to the data for the packing fraction at 
jamming for a > 0 . 8 . (b) Rc plotted versus the normalized 
difference Apj/(pj), where Apj = tpj — pj. The master 
curve (solid line) obeys logj^Q Ac {Apj/ [T^- In (n) 
and (b), we considered binary systems with diameter ratios 
a > 0.8 (crosses) and ternary systems with diameter ratio 
pairs gbJga = 0.95 and gcIga = 0.9, gb/ga = 0.9 and 
gcIga ~ 0.88, and gbIga = 0.95 and gc/ga ~ 0.5 (squares, 
circles, and triangles, respectively). In the absence of demix¬ 
ing, Ac versus the jammed packing fraction for ternary sys¬ 
tems is quantitatively similar to that for binary systems. 
However, ternary systems that demix deviate from the master 
curve. 

creases the particle size gradient in the original binary 
system without changing the mechanism that drives crys¬ 
tallization [ 1 ^ [ 23 . The deviations from the master curve 
(solid lines) in Fig. [S] indicate demixing in the ternary 
system with diameter ratio pairs gb/ga = 0.95 and 
xc/xa = 0.5, and in this case higher packing fractions 
are obtained than for binary systems. 

In Fig. [HI we illustrate the diameter ratio variation for 
15180 combinations of three elements from a set of 46 
potential BMG-forming elements. The atomic sizes of 
the elements are given by their metallic radii Bin, 
which are shown in the inset of Fig. [0] We define the 
component types so that they satisfy gc gb ^ ga, 
and thus the data occurs in left corner of the plot. 

In previous studies, we predicted that the optimal bi¬ 
nary hard-sphere glass formers occur in the diameter ra¬ 
tio range 0.73 < a < 0.82, where a is sufficiently small 
to prevent ordering, but not too small to cause demixing. 
We can estimate the optimal glass-forming regime in the 
xb/xa and gc/xa plane for ternary systems using the 
following arguments. First, two of the boundaries {ab 
and cd) can be obtained directly from the results for bi¬ 
nary systems. Because adding a third intermediate-sized 
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component can prevent demixing of the original two com¬ 
ponents, we expect that the lower bound for the diam¬ 
eter ratio in ternary BMG-forming systems to be much 
smaller than 0.73. We propose that 0.73^ is the lower 
bound for the diameter ratio for ternary systems. In this 
case acjoA = 0.73^, gbIcfa = 0.73, and acl<JB = 0.73 
(point e in Fig. 15]), and thus all binary combinations are 
above the lower bound of the good GFA regime. 

We predict that good BMG-forming alloys will oc¬ 
cur within the polygon defined by lines connecting the 
points (a)-(e) in the crs/o-yi and acjoA parameter space. 
33 ternary alloy systems have been observed exp erimen- 
tally in amorphous states (filled circles) @,[11 Hi, all 
of which fall in the good glass-forming regime predicted 
by hard-sphere systems. For the experimentally ob¬ 
served ternary BMGs, the diameter ratio for the small¬ 
est versus the largest particle acjoA satisfies a < 0.8, 
which is below the demixing limit for binary systems. 
Therefore, the experimentally observed ternary BMGs 
have better GFA than the best binary glass-forming al¬ 
loys. In addition, the experimentally observed BMGs 
tend to be positioned away from the boundaries ab and 
cd. As ternary systems approach ab (cd), obI<Ja 1 
[pcI’^A 1), which causes them to behave as binary 
systems and reduces their glass-forming ability. The ex¬ 
perimentally observed ternary BMGs cluster roughly into 
three groups: (i) (Zr,Hf,Sn,Mg)-(Al,Ti,Nb)-(Gu,Ni,Co), 
(ii) (Y,Ln)-Al-(Cu,Ni,Go), and (iii) (Au,Pd,Pt)-(Cu,Ni)- 
(Si,P). (See Table |T] in Appendix [Ml 


4. CONCLUSION 

We performed event-driven molecular dynamics simu¬ 
lations of ternary hard-sphere systems over the full range 
of compositions and diameter ratios to identify the op¬ 
timal glass-forming parameter regime. We identify two 
mechanisms for optimizing the glass-forming ability in 
ternary systems. First, if the sizes of the three compo¬ 
nents are similar, i.e. less than 10% deviation in the 
diameters, the ternary system behaves effectively as a 
binary system containing only the largest and smallest 
particles. Second, if the diameter ratio of the small¬ 
est to the largest particle is below the demixing thresh¬ 
old ac/cFA ^ 0.8, adding the third component B with 
A < (^bI'aA < 1 will dramatically enhance the glass¬ 
forming ability above that for binary systems. We show 
that all experimentally observed ternary BMGs to date 
possess atomic species for which the diameter ratio of the 
smallest to the largest satisfies acjcfA < 0.8. Thus, an 
efficient strategy to design BMGs with good GFA is to 
maintain large atomic size differences and prevent demix¬ 
ing by introducing three or more atomic components. 

We recognize that the inter-atomic potentials describ¬ 
ing BMGs are far more complicated than the pairwise 
additive hard-sphere potential that we employed. For ex¬ 
ample, the apparent distance between the repulsive cores 
of two elements can be shorter than the mean core size of 



ac/aA 

FIG. 6: (Color online) Scatter plot of the diameter ratios 
ns/CTA versus crc/crA with ac < us < cta for all 15180 com¬ 
binations of three elements chosen from 46 possible BMG- 
forming elements (dots). 33 of the combinations (filled circles) 
have been shown experimentally to form BMGs, and occur 
in roughly three main clusters: i) (Zr,Hf,Sn,Mg)-(Al,Ti,Nb)- 
(Cu,Ni,Co), ii) (Y,Ln)-Al-(Cu,Ni,Co), and iii) (Au,Pd,Pt)- 
(Cu,Ni)-(Si,P). The predicted BMG-forming alloys are lo¬ 
cated within the polygon bounded by the solid lines. The 
inset gives the atomic radii in pm of the 46 potential BMG- 
forming elements ordered from smallest to largest [sl. fl^. fl^ . 
The symbols of the elements alternate from top to bottom on 
the horizontal axis and the lanthanide elements are shown as 
filled circles. 

the two elements [l^. Also, in a more exact treatment, 
Friedel oscillations originating from perturbations to the 
electron density should be included since they are known 
to change the stability of various crystalline lattices [s^ ■ 
Despite these caveats and others, we show that the hard- 
sphere model is able to semi-quantitatively predict the 
regime of optimal glass-forming ability in experimentally 
observed ternary BMGs. In the near future, we will con¬ 
sider how non-additivity of the particle diameters, at¬ 
tractive interactions and barriers in the pairwise poten¬ 
tial, and multi-body interactions affect crystal and glass 
formation and modify the hard-sphere predictions [3l|. 
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Appendix A: Glass-forming ability of experimentally 
observed BMGs 

In this appendix, we provide a table of the critical cool¬ 
ing rates TZc in units of K/s for experimentally observed 
binary and ternary BMGs. See Table HI In the first col¬ 
umn, we list each class of binary and ternary BMGs ac¬ 
cording to the atom types that are present. Atom types 
with similar sizes and properties are grouped together in 
parentheses. In the second column, we provide examples 
of specific alloys within each BMG class. The third col¬ 
umn gives diameter ratios: ubIcfa with as < cfa for bi¬ 
nary systems and (JbIo a and an la a with ac < as < aA 
for ternary systems jl, |T3, [H, [1^ . 

Appendix B: Measurement of Critical Compression 
Rate Rc in Simulations 

In this appendix, we describe the measurement of the 
critical compression rate Rc in the molecular dynam¬ 
ics simulations of hard-spheres. We performed 96 inde¬ 
pendent runs to generate jammed configurations at each 
compression rate R. We find that the distribution P{Qe) 
of the global bond orientational order parameter becomes 
bimodal as R ^ Rc with peaks that correspond to amor¬ 
phous and partially crystalline configurations. In Fig. [3 
we show that both the mean and median Qg possess a 


sigmoidal shape on a logarithmic scale in R. We define 
Rc as the critical compression rate at which the mean 
and median Qq intersect. A configuration is determined 
to be crystalline (amorphous) if Qe > Qc {Qe < Qc), 
where Qc is the value of Qe at which the mean and me¬ 
dian Qq intersect. Since the distribution P{Qq) is bi¬ 
modal at i? w R^^ one can also fit P{Qq) by the sum of 
two Gaussian distributions and identify Rc as the rate at 
which the two Gaussian contributions have equal area. 
This method can be more time efficient since it avoids 
measurements at i? -C i?c- 



FIG. 7: (Color online) The mean and median Qq as a function 
of compression rate R for ternary hard spheres with diameter 
ratios acjaA = 0.88 and obIoa = 0.9 and compositions 
xa = xb = xc = 1/3. The critical compression rate Rc 
and bond orientational order parameter Qc are defined by 
the intersection of the mean and median Qq. The solid lines 
are fits to a logistic function. 
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TABLE I: Glass-forming ability (characterized by the critical cooling rate TZc in units of K/s) and 
atomic diameter ratio“ for experimentally observed binary and ternary BMGs (sl. [T^ [^. Elements 
with similar sizes and properties are grouped together using parentheses (see also Fig. [I. The 
lanthanide elements*’ are indicated by Ln. 


Binary system 

Alloy 

77c (K/s) 

Diameter ratio (Js/cta 

Fe-B 

FegiBg 

2.6 X 10' 

0.62 

(Au,Pd)-Si 

AusoSigo 

1.0 X 10® 

0.71 


PdgsSis 

5.0 X 10'^ 

0.74 


PdsgSiis 

1.8 X 10® 



PdrsSigs 

1.0 X 10® 


Ti-Be 

Ties Bear 

6.3 X 10® 

0.76 

Zr-Be 

ZresBess 

1.0 X 10^ 

0.7 

Zr-Cu 

ZrsoCuso 

250 

0.8 

Nb-Ni 

Nb4oNi6o 

1400 

0.85 

Ternary system 

Alloy'’ 

77c (K/s) 

Diameter ratios ctc/ua, ctb/ua 

Au-Si-Ge 

Au77,8Si8.4Gei3.8 

3.0 X 10® 

0.71, 0.79 

(Au,Pd,Pt)-(Cu,Ni)-(Si,P) 

Pd4oNi4oP20 

0.167 

0.73, 0.91 


Pd77CU6Sil7 

125 

0.74, 0.93 


Pd7g.5CU4Sil6.5 

500 



Pd77.5CU6Sil6.5 

100 


(Y,Ca,Ln)-Mg-(Cu,Ni) 

Ndi5Mg7oNii5 

178.2 

0.69, 0.88 


NdisMgesNigo 

30 



NdioMg75Nii5 

46.1 



Nd5Mg77Nii8 

49000 



NdsMgggNis 

53000 



NdioMg8oNiio 

1251.4 



YioMgesGugs 

50 

0.71, 0.89 


GdioMges Cu25 

1 

0.71, 0.89 

(Y,Ln)-Al-(Gu,Ni,Co) 

La55Al25Ni2o 

67.5 

0.66, 0.76 


La55Al25Cu20 

72.3 

0.68, 0.76 


La66All4Cu20 

37.5 


(Zr,Hf,Sn,Mg)-(Al,Ti,Nb)-(Cu,Ni,Co) 

ZreeAlsNige 

66.6 

0.78, 0.89 

(Zn,Al,Ag)-Mg-Ga 

ZnggMgisGaes 

20 

0.69, 0.81 


^ Atomic radii are obtained from Refs. [l^. [T^. which determine the atomic sizes using the first 
peak of the radial distribution function of amorphous liquid alloys or half of the spacing between 
atoms in metallic solids. 


^ Ln refers to the series of fifteen metallic elements (La, Ge, Pr, Nd, Pm, Sm, Eu, Gd, Tb, Dy, Ho, 
Er, Tm, Yb, and Lu) with atomic numbers 57-71. Together with two more chemically similar 
elements Sc and Y, these seventeen elements are collectively known as the rare earth elements 
and are typically the largest sized component in BMGs. 

For each system, we only list alloys for which the critical cooling rate TZc has been reported. Other 
alloys such as Ga-Mg-Cu, Hf-Al-Cu, and Y-Al-Go, are also BMG formers, but with unreported 
values of 77c. 
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